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Abstract: This paper offers a geometric framework for modeling port-Hamiltonian systems on 
discrete manifolds. The simplicial Dirac structure, capturing the topological laws of the system, 
is defined in terms of primal and dual cochains related by the coboundary operators. This finite- 
dimensional Dirac structure, as discrete analogue of the canonical Stokcs-Dirac structure, allows 
for the formulation of finite-dimensional port-Hamiltonian systems that emulate the behaviour 
of the open distributed-parameter systems with Hamiltonian dynamics. 
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1. INTRODUCTION 

A large class of open distributed-parameter Hamiltonian 
systems can be defined with respect to the Stokes-Dirac 
structure van der Schaft and Maschke [2002] . This infinite- 
dimensional Dirac structiirc provides a theoretical account 
that permits the inclusion of varying boundary variables 
in the boundary problem for partial differential equations. 
From an interconnection and control viewpoint, such a 
treatment of boundary conditions is essential for the incor- 
poration of energy exchange through the boundary, since 
in many applications the interconnection with the envi- 
ronment takes place precisely through the boundary. For 
numerical integration, simulation and control synthesis, it 
is of paramount interest to have finite approximations that 
can be interconnected to one another. 

Most of the numerical algorithms emanating from the field 
of numerical analysis and scientific computing, however, 
fail to capture the intrinsic system structures and proper- 
ties, such as symplecticity, conservation of momenta and 
energy, as well as differential gauge symmetry. Further- 
more, some important results, including the Stokes theo- 
rem, fail to apply numerically and thus lead to spurious 
results. 

Recently in Seslija et al. [2011a], we have suggested a dis- 
crete exterior geometry approach to structure-preserving 
discretization of distributed-parameter port-Hamiltonian 
systems. The spatial domain in the contimious theory 
represented by a finite-dimensional smooth manifold is 
replaced by a homological manifold-like simplicial complex 
and its circumcentric dual. The smooth differential forms, 
in discrete setting, are mirrored by cochains on the primal 
and dual complexes, while the discrete exterior derivative 



is defined to be the coboundary operator. A discrete ana- 
logue of the Stokes-Dirac structure is a so-called simplicial 

Dirac structiu-c defined over a space of primal and dual 
discrete differential forms. 

In this paper, we address the issue of matrix repre- 
sc!ntations of simplicial Dirac structures by representing 
cochains by their coefficient vectors. In this manner, all lin- 
ear operator from the continuous world can be represented 
by matrices, including the Hodge star, the coboundary 
and a trace operator. First, we recall the definition of the 
Stokes-Dirac structme and port-Hamiltonian systems. In 
the third section, we define some essential concepts from 
discrete exterior calculus as developed in Desbrun et al. 
[2002, 2003], Hirani [2003]. In order to allow the inclusion 
of nonzero boundary conditions on the dual cell complex, 
in Seslija et al. [2011a] we have adapted a definition of the 
dual boundary operator that leads to a discrete analogue of 
the integration by parts formula, which is a crucial ingredi- 
ent in establishing simplicial Dirac structures on a primal 
simplicial complex and its circumcentric dual. Finally, we 
demonstrate how these simplicial Dirac structures relate to 
spatially discrctized wave equation on a bounded domain 
and the telegraph equations on a segment. 



2. THE STOKES-DIRAC STRUCTURE AND 
PORT-HAMILTONIAN SYSTEMS 

The Stokes-Dirac structure is an infinite-dimensional Dirac 
structure that provides a foundation for port-Hamiltonian 
formulation of a class of distributed-parameter systems 
with boundary energy flow van der Schaft and Maschke 
[2002]. 



Throughout this paper, let M be an oriented n-dimensional 
smooth manifold with a smooth (n — l)-dimensional 
boundary dM endowed with the induced orientation, rep- 
resenting the space of spatial variables. By f2'^(M), k = 
0, 1, . . . , n, denote the space of exterior fc-forms on M, and 
by rj'^(9M), A; = 0, 1, . . . ,n — 1, the space of fc-forms on 
dM. 

For any pair p, q of positive integers satisfying p-\-q = n+1, 
define the flow and effort linear spaces by 

Fp^q =QF{M) X n^iM) X Q"-f(9M) 

£p^q =n"-P{M) X f7"-«(M) X n"-'?(9M) . 

The bilinear form on the product space J-p^q x £p^q is 

(({/p T fq T fb ^ ^p' ^q' ^fe)' {fp ^ fq ^ fb ^ ^p' ^g' ^b))) 
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Theorem 1. Given linear spaces Tp^q and Sp^q, and the 
bilinear form, ((,)), define the following linear subspace V 
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where d is i/ie exterior derivative and {qm stands for a 
trace on the boundary dM. Then V = T>^ , that is, V is a 
Dirac structure. 



Consider a Hamihonian density H : VLP{M) x ri«(M) -J> 
n"(M) resulting with the Hamihonian H = /^'H € M. 

Setting the flows fp = fq = ^^'^ efforts 

Cp = SpH, Cq = SqH, where (SpH.SqH) e f7"-P(M) x 
f2"~'(M) are the variational derivatives of H at {ap,aq), 
the distributed-parameter port-Hamiltonian system is de- 
fined by the relation 



dap 



da. 
'~di 



^,fb,6pH,6qH,ebj gV, t 



Since ^ = eb A fb, the system is lossless. 



3. BASICS OF DISCRETE EXTERIOR CALCULUS 

In the discrete setting, the smooth manifold M is re- 
placed by an oriented manifold-like simplicial complex. 
An n-dimensional simplicial complex ii' is a simplicial 
triangulation of an n-dimensional polytope \K\ with an 
(n — l)-dimensional boundary. Familiar examples of such 
a discrete manifold arc meshes of triangles embedded in 
and tetrahedra obtained by tetrahedrization of a 3- 
dimensional manifold. 

3.1 Chains and cochains 

The discrete analogue of a smooth /c-form on the manifold 
M is a A;-cochain on the simplicial complex K. A k- 
chain is a formal sum of fc-simplices of K such that 



its value on a simplex changes sign when the simplex 
orientation is reversed. The free Abelian group generated 
by a basis consisting of oriented fc-simplices with real- 
valued coefficients is Ck{K;M.). The space Ck{K;R) is a 
vector space with dimension equal to the number of fc- 
simplices in K, which is denoted by Nk- The space of k- 
cochains is the vector space dual of Ck{K;R) denoted by 
C'^{K;'R) or ^l'^{K), as a reminder that this is the space 
of discrete fc-forms. The space ^^iK) is the space of real- 
valued linear functionals on the vector space Ci.{K;R). 

The discrete exterior derivative d'^ : C'^{K) C''~^^{K) 
is defined by duality to the boundary operator dk+i ■ 
Ck+i{K;Z) — > Ck{K;Z), with respect to the natural 
pairing between discrete forms and chains. For a discrete 
form a e {K) and a chain Ck+i € Ck+i {K; Z) we define 
d^^ by 

(d a,Cfc+i) = (a, (d )'^Cfc+i) = {a,dk+iCk+i) , 

where the boundary operator d^+i is the incidence matrix 
from the space of A;-|-l-simplices to the space of fc-simplices 
and is represented by a sparse N^+i x Nk matrix containing 
only or ±1 elements Desbrun et al. [2008]. The important 
property of the boundary operator is dk o dk+i = 0. 
The exterior derivative as the coboundary operator also 
satisfies d^'^^ o d*^ = 0, what is a discrete analogue of the 
vector calculus identities curlograd = and divocurl = 0. 



(2) 3.2 Dual cell complex 



An essential ingredient of discrete exterior calculus is the 
dual complex of a manifold-like simplicial complex. Given 
a simplicial well-centered complex K , we define its interior 
dual cell complex -kiK (block complex in terminology of 
algebraic topology Munkres [1984]) as a circumcentric 
dual restricted to An important property of the the 
Voronoi duality is that primal and dual cells are orthogonal 
to each other. The boundary dual cell complex -k^K is 
a dual to dK. The dual cell complex -kK is defined as 
■kK = kiK X k\jK. A dual mesh -k^K is a dual to K in sense 
of a graph dual, and the dual of the boundary is equal to 
the boundary of the dual, that is d{-kK) = -k(dK) ~ -k^K. 
Because of duality, there is a one-to-one correspondence 
between fc-simplices of K and interior (n — fc)-cells of -kK. 
Likewise, to every fc-simplex on dK there is a uniquely 
associated (n — 1 — fc)-cell on d{-kK). Fig. 1 illustrates the 
duality on a flat 2-dimensional simplicial complex. 

In order to properly account for the behaviours on the 
boundary, we need to adapt the definition of the boundary 
dual operator as presented in Hirani [2003] , Desbrun et al. 
[2003]. In Scslija et al. [2011a] we have proposed the 
following definition. 

Definition 1. The dual boundary operator dk ■ Cfe(*iii';Z) 
— Ck-i{*K;Z,) is a homomorphism defined by its action 
on a dual cell a'^ = ★icr""'^ = *i[t;o, • • • ) Vn-k], 

dka'' = d M^k [VO,---, Vn-k] = di,k *i [l^O, • • • , Vn-k] 

+ dh,kM [vq,. ■ . ,Vn-k], 

where 

di,k *i [vo, Vn-k] = X] *i(s^„-fc+i(T"'~''"'"^) 

dh,k *i [vo, Vn-k] = *b {s^n-ka'^~'') . 




Fig. 1. A 2-diniensional simplicial complex K and its cir- 
cumcentric dual cell complex -kK indicated by dashed 
lines. The boundary of -kK is the dual of the boundary 
of i^. 

The boundary of the dual cell complex as defined in Hirani 
[2003] is equal to di. The dual boundary operator dh 
extends the definition from Hirani [2003] in such a manner 
that the boundary of the extended dual cell complex -kK 
is the geometric boundary. The dual exterior derivatives 
df-i : C^-^{-kiK) C^{-k;K) and d^-^ : C^-^(:k^K) ^ 
C''{-kiK) are defined by duality to the dual boundary 
operators di^k and dh,k, respectively. 

3.3 Discrete wedge and Hodge operator 

There exists a natural pairing, via the so-called primal- 
dual wedge product, between a primal fc-cochain and a 
dual {n — fc)-cochain. The resulting discrete form is the 

volume form. Let a'' € ^ai^) and /3"-'= e ny''{*iK). 
We define the discrete primal-dual wedge product A : 

n'^^iK) X n",-'{*iK) ^ n^{Vk{K)) by a V,.) = 

(a'=,cr'=)(/3"-'=,*iCT'=) = (-l)'''("-'=)(/3"-'= Aa'=,K,fc), where 
V^k is the n-dimensional support volume obtained by 
taking the convex hull of the simplex a'^ and its dual -ka^. 

The proposed definition of the dual boundary operator 
ensures the validity of the evaluation by parts relation 
that parallels the integration by parts formula for smooth 
differential forms. 

Proposition 2. Let K he an oriented well- centered sim- 
plicial complex. Given a primal [k — l)-form a and an 
internal dual (n — k)-discrete form /3i e Vl"^^^ {k-^K) and a 
dual boundary form /3b G ^^~^{'^hK), then 

(d'^+^a A A, if) + (-f VA(dr-'=/3i + d^-'=/3b), K) 

= {tY^-^ahP^,dK) , 



where tr 



fe-i 



is the trace operator that isolates the com- 



ponents of the primal {k — l)-cochain associated with the 
boundary dK . 



From definition of the dual boundary operator follows 
that d['-'= = (-f)'=(d'=-i)^ and d'^^'' = {-l)''-\tr''-^y 
(confer to Seslija et al. [20f lb]). 

The support volumes of a simplex and its dual cell are the 
same, which suggests that there is a natural identification 
between primal fc-cochains and dual (n — fc)-cochains. In 
the exterior calculus for smooth manifolds, the Hodge star, 
denoted is an isomorphism between the space of Re- 



forms and (n — A:)-forms. The discrete Hodge star is a 

map *k ■ ^di^) ~^ ^d~''i*i^) defined by its value over 
simplices and their duals. The Hodge star *fc is a diagonal 
Nk X Nk matrix with the entry corresponding to a simplex 
a'' being |cr'^|/| ★ cr'^|. 

Another possibility for the construction of the Hodge oper- 
ator us to use Whitney forms. The Whitney map is an in- 
terpolation scheme for cochains. It maps discrete forms to 
square integrable forms that are piecewise smooth on each 
simplex. The Whitney maps are built from barycentric 
coordinate functions and the resulting matrix is sparse but 
in general not diagonal Bossavit [1998], Hiptmair [2002]. 

4. SIMPLICIAL DIRAC STRUCTURES 

In this section, we introduce Dirac structures defined in 
terms of primal and duals cochains on the underlying 
discrete manifold. We call these Dirac structures simplicial 
Dirac structures. 

In the discrete setting, the smooth manifold M is replaced 
by an n-dimensional well-centered oriented manifold-like 
simplicial complex K. The flow and the effort spaces will 
be the spaces of complementary primal and dual forms. 
The elements of these two spaces are paired via the discrete 
primal-dual wedge product. Let 



and 



^T^iK) X n'l-^k.K) X n'l-^dii^K)) . 



The primal-dual wedge product ensures a bijective relation 
between the primal and dual forms, between the flows and 
efforts. A natural discrete mirror of the bilinear form (1) 
is a symmetric pairing on the product space iF^^^ x S^ ^ 
defined by 



{{{ fp ^ fq T fb ^ 




(f2 r2 f2 

\J p 1 J q 1 J b 1 



= {el A fl + e\h fl + elA fl + a f^K) 
+ {elhfl + elhfldK). 



(3) 



A discrete analogue of the Stokes-Dirac structure is the 
finite-dimensional Dirac structure constructed in the fol- 
lowing theorem Seslija et al. [2011a]. 

Theorem 3. Given linear spaces ^ and £p^q, and the 
bilinear form {{,))d- The linear subspace Vd C ^ x E^^^ 
defined by 

T^d = {(/pi fqi fbi Gp, iq, Pft) G X £p q\ 

(-l)Mr''^ ^ /^"--j 







+ (-ir 



dri— 
b 





fb = (-l)^'tr"-Pep} , 



with r — pq - 
pairing {{,))d 



(4) 

1, is a Dirac structure with respect to the 



Note that since 
(_l)"-p(tr"-P)T 



^ (-l)9(d"-P)^ 
the structure (4) is in fact a Poisson 



and d[^' 



structure on the state space n^{k^K) x fl'^{K). 

The other discrete analogue of the Stokes-Dirac structure 
is defined on the spaces 



-\dK). 

A natural discrete mirror of (1) in this case is a symmetric 
pairing defined by 

((( /p I /(J 1 /b J Cp) Cq; Cfc)) {fp J fqJ fb' ^p' ^g' ^b))) d 



= (ep A fl + ej A + A fl + A Z^, Jf) 
+ (e^A/2 + egA/^ajr). 
Theorem 4. T/ie linear space Vd defined by 

i'd = {{.fpj fq- .fhTep,eq,eb) G X fp ql 



+ 







fb, (5) 



(-l)P9+ld"-«^ 

VC"^ V 

et = (-l)Ptr"-9e,} 
is a Dirac structure with respect to the bilinear pairing 

Ud- 

In the following section, the simplicial Dirac structures 

(4) and (5) will be used as terminus a quo for the geo- 
metric formulation of spatially discrete port-Hamiltonian 
systems. 

5. PORT-HAMILTONIAN SYSTEMS ON A 
SIMPLICIAL COMPLEX 



An alternative formulation of a spatially discrete port- 
Hamiltonian system is given in terms of the simplicial 
Dirac structure (5). We start with the Hamiltonian func- 
tion (ap,dg) H-^ 'H(ap,dg), where Up G ^^d^K) and 
dq e Q!j^{i^iK). In a similar manner as in deriving (7), 
we introduce the port-Hamiltonian system 

/ dap 



\ dt 



(- 



1 







Bb = (-l)Ptr"- 





fdn\ 




dap 




&H 




\daqj 


on 




d&q ' 





+ 



fb , 



(8) 



In contrast to (7), in the case of the formulation (8), 
the boundary flows can be considered to be freely 
chosen, while the boundary efforts are determined by 
the dynamics. The free boundary variables are always 
defined on the boundary of the dual cell complex. 

6. PHYSICAL EXAMPLES 

In this section we consider the discrete wave equation 
on a 2-dimensional simplicial complex and the telegraph 

equations on a segment. 

6.1 Two-dimensional wave equation 



Let a function H : ^^i^iK) x n^K) R stand for the 
Hamiltonian (dp, a,) H{ap,aq), with dp G Q^{*iK) 
and Uq e Q,^{K). A time derivative of H along an arbitrary 
trajectory t {ap{t) , aq{t)) € 0,li*iK) x O^(iV'), i e R, is 



—n{ap,aq) 



da„ 



A 



Ik 



dn 

da„ 



A 



daq 

1h 



.K). (6) 



The relations between the simplicial-Dirac structure (4) 
and time derivatives of the variables are: fp = —^§f, 



da 



g^r, while the efforts are: ep 



dH 

d&r, 



dH 

daa ' 



This allows us to define time-continuous port-Hamiltonian 
system on a simplicial complex K (and its dual -kK) by 




_f (-i)'-dr 

d"-f 



f—\ 

dap 
dH 
\daqj 



+ (-1)' 



(7) 



(-Iftr" 



, dU 

d&r, 



where r = pq + 1. 
It immediately follows that 



d u 



{eb A fb, dK), enunci- 



ating a fundamental property of the system: the increase 
in the energy on the domain \K\ is equal to the power 
supplied to the system through the boundary dK and 
d{-kK). The boundary efforts Cb are the boundary control 
input and fb are the outputs. 

Remark 1. Introducing a linear negative feedback control 
as e-b = (— where *b is the Hodge star 
on the boundary dK , leads to passivization of the lossless 
port-Hamiltonian system, i.e., ^H < —{fb A *hfb,dK) < 
0. Furthermore, if the Hamiltonian is a fCoc, function with 
a strict minimum that is a stationary set for the system 
(7), the equilibrium is asymptotically stable. 



Let us consider the simplicial Dirac structure underpining 
the discretized two-dimensional wave equation. The nor- 
malized wave equation is given by 

where (p'^ is a smooth 0-form on a compact surface M cM.^ 
with a closed boundary, and A is the Laplace operator. 
(Throughout, the superscript c designates the continuous 
quantities.) This equation, together with nonzero energy 
flow, can be formulated as a port-Hamiltonian system with 
boundary port variables Golo [2004], Talasila [2002]. 

The energy variables of the discretized system are chosen 
as follows: the kinetic momentum is a dual 2-form whose 
time derivative is set to be fp, the elastic strain is a 
primal 1-form with time derivative corresponding to fq, 
the coenergy variables are a primal 0-form Cp and a dual 
1-form Cq. Such a formulation of the discrete wave equation 
is consonant with the simplicial Dirac structure (4) for the 
case when p = n = 2 and q — 1, and is given by 

-dl 






dO 







fb 



tr^e 



p • 



Example. Consider a simplicial complex pictorially given 
by Fig. 2. The primal and dual 2-faces have counter- 
clockwise orientations. The matrix representation of the 
incidence operator di , from the primal edges to the primal 
vertices, is 



V2 
V3 



[V0,Vl\ 

-1 

1 






[Vl,V2] 



-1 
1 




[V2,Vo\ 




-1 
1 



[Vl,V3\ 



-1 


1 



[V3,V2\ 





1 

-1 




«0 



V2 



V2n-2 



7 V3 



Fig. 2. A simplicial complex K consists of two triangles. 
The dual edges introduced by subdivision are shown 
dotted. 



while the discrete exterior derivative from the vertices to 
the edges is the transpose of the incidence operator, i.e., 
— dj. The dual exterior derivative is d?^ = — (d°) , 
while the matrix representation of the dh.2 operator is 

★iWO ★iUl -*iV2 *iV3 

[v2,vi] 10 

[VUV3] 10 

[V3,V4] 1 

[V4,V2] 10 

The trace operator is tr*^ = i'^hV — dh,2- 
It is trivial to show 

(dOgp A e„ K) + (ep A (d^e, + diet), K) 
= eb[v2,vi]fb{vo) + eb[vi,vz]h{vi) 
+ eb[v3,Vi\h[v3) + eb[vi,V2]fb{v2) , 



(9) 



what confirms that the boundary terms genuinely live on 
the boundary of \K\. 



6.2 Telegraph equations 

We consider an ideal lossless transmission line on a 1- 
dimensional simplicial complex. The energy variables are 
the charge density q £ 17j(A'), and the flux density (f> S 
rij(*if), hence p = q = 1. The Hamiltonian representing 
the total energy stored in the transmission line with 
distributed capacitance C and distributed inductance L 
is 



n 



1 



1 



(10) 



where *o and *i are the discrete diagonal Hodge operators 
that relate the appropriate cochains according to the 
following schematic diagram 



1*0 



1*1 



The co-energy variables are: e 
representing voltages and Cg 



-p " 

dU 
d4> 



d-H _ ^ 9 _ 



= V 



I currents. 



Selecting /p = -|2 and fq = 
Hamiltonian formulation of the telegraph equations 



^ leads to the port- 



V2n-1 V2t, 



Fig. 3. The primal 1-dimensional simplicaial complex K. 
By construction, the nodes vq and V2n are added to 
the boundary to insure that d{-kK) = -k{dK). 




d" 
d? 



/ 



° L 



(11) 



where fh are 
currents. 



the input voltages and Cb are the output 



In the case we wanted to have the electrical currents as the 
inputs, the charge and the flux densities would be defined 
on the dual mesh and the primal mesh, respectively. 
Instead of the port-Hamiltonian system in the form (11), 
the discretized telegraph equations would be in the form 
(7). The charge density is defined on the dual cell complex 
as g G Vl\{*iK) and the discrete flux density is </) G ^\{K). 
The finite-dimensional port-Hamiltonian system is of the 
form 



/ 



h - -tr" 




(12) 



where Ch are the input currents and fb are the output 
voltages. 

The exterior derivative d° : ^°^{K) ^\{K) is the 
transpose of the incidence matrix of the primal mesh. The 
discrete derivative d? : 0^(*iX) — Vl]i{-k{K) in the matrix 
notation is the incidence matrix of the primal mesh. Thus, 
we have 



0^T 



tr" = {dlY = 



/-I 1 
' 0-11 



V 

-1 •• 
•• 



0\ 


-1 1/ 



(13) 




1 



Remark 2. The discrete analogue of the Stokes-Dirac 
structure obtained in Golo [2004] finite- dimensional 
Dirac structure, but not a Poisson structure. The impli- 
cation of this on the physical realization is that the trans- 
mission line in the finite- dimensional case is not only com- 
posed of inductors and capacitors but also of transformers. 

The physical realizations of the port-Hamiltonian systems 
(11) and (12) are given on Fig. 4 and Fig. 5, respectively. 
Stabilization of either of those systems is easily achieved by 
terminating boundary ports with resistive elements, what 



— nm^^—r-Jim^^—r-rm^'''-^ 

n, —I— rr. —I— r7„ — L 



T T 



Fig. 4. The finite-dimensional approximation of the lossless 
transmission line when the inputs are voltages and the 
outputs are currents. The inductances Li,...,L„+i 
are the values that the discrete distributed inductance 
L takes on the simplices [vg, wi], . . . , [w2n-i, i'2n\y the 
capacitances Ci , . . . C„ are the values C takes on 
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Fig. 5. The finite-dimensional approximation of the lossless 
transmission line when the inputs are currents and 
the outputs are voltages. The inductances are: Li — 

hvoM^' = ^(N,V2]), L2 = L-" = L{[v2,Vi]), 

■ ■ •' ^" " /[f2„-2.i'2„] = L([v2n-2,V2n]); the valucs 
of capacitors are: Ci — /j.^ .^j — C{[vo,vi]), C2 = 

■ ■ '^"+1 = /[fi2„-l,fi2„] = C{[v2n-l,V2n])- 

is a practical application of the passivization explained in 
Remark 1. 

The accuracy of the proposed method is 1/n (see Seslija 
et al. [2011b]). 
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